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ABSTRACT 

Terrestrial planets in the solar system, such as the Earth, are oxygen-rich, with silicates and iron 
being the most common minerals in their interiors. However, the true chemical diversity of rocky 
planets orbiting other stars is yet unknown. Mass and radius measurements are used to constrain 
the interior compositions of super-Earths (cxoplancts with masses of 1 - 10 M^), and are typically 
interpreted with planetary interior models that assume Earth-centric oxygen-rich compositions. Using 
such models, the super-Earth 55 Cancri e (mass 8 Af^, radius 2 R^) has been suggested to bear an 
interior composition consisting of Fc, silicates, and an envelope (> 10% by mass) of super-critical 
water. We report that the mass and radius of 55 Cancri e can also be explained by a carbon-rich solid 
interior made of Fe, C, SiC, and/or silicates and without a volatile envelope. While the data allow 
Fe mass fractions of up to 40%, a wide range of C, SiC and/or sihcate mass fractions are possible. A 
carbon-rich 55 Cancri e is also plausible if its protoplanetary disk bore the same composition as its 
host star, which has been reported to be carbon-rich. However, more precise estimates of the stellar 
elemental abundances and observations of the planetary atmosphere are required to further constrain 
its interior composition. The possibility of a C-rich interior in 55 Cancri e opens a new regime of 
geochemistry and geophysics in extraterrestrial rocky planets, compared to terrestrial planets in the 
solar system. 

Subject headings: planetary systems — planets and satellites: general — planets and satellites: indi- 
vidual (55 Cancri e) 



1. INTRODUCTION 

The recently discovered transiting super-Earth 55 Can- 
cri e orbiting a nearby (naked eye) G8V star every 18 
hours is one of the most favorable super-Earths for de- 
tailed characterization. The planet was originally discov- 
ered in a five-planet system by the radial velocity method 
which gave its minimum mass (McArthur et al. 2004; 
Fischer ct al. 2008). However, the recent discovery of its 
transits across the host star (Winn et al. 2011; Demory 
et al. 2011) has made possible accurate measurements 
of its mass (8.37 ± 0.38 M®; Endl et al. 2012) and ra- 
dius (2.04 ± 0.15 Re; Gillon ct al. 2012), and, hence, 
constraints on its interior composition. 

The interior of 55 Cancri e, inferred from its mass and 
radius, has been suggested to bear an interior composi- 
tion made of iron, silicates, and an envelope (> 10% by 
mass) of super-critical water (Winn et al. 2011, Demory 
et al. 2011; Gillon et al. 2012). A water envelope is 
required to explain the mass and radius with an oxygen- 
rich interior composition. Given its mass, even a purely 
silicate interior with no iron content could only barely 
explain the observed radius. Consequently, a massive 
water envelope (> 10% by mass) over a rocky Earth-like 
interior has been required in the above studies to explain 
the radius of this highly irradiated super-Earth. More- 
over, given the extreme stellar irradiation (T ^ 2500 K; 
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Demory et al. 2012), the water has been suggested to be 
in a supercritical state in the planetary envelope. 

The silicate-rich and water-rich interior of 55 Cancri e 
as suggested by previous studies is based on the assump- 
tion of an oxygen-rich interior, motivated by terrestrial 
planets in the solar system. However, several recent stud- 
ies have suggested the possibility of carbon-rich extraso- 
lar planets, both theoretically (Lodders 2004; Kuchner 
& Seager 2005; Bond et al. 2010; Madhusudhan et al. 
2011b; Obcrg ct al. 2011; Mousis ct al. 2012) and ob- 
scrvationally (Madhusudhan ct al. 2011; Madhusudhan 
2012). Motivated by these studies, in the present work 
we investigate if the observed mass and radius of 55 Can- 
cri e may be explained by a carbon-rich interior. 

The possibility of a carbon-rich interior in 55 Cancri e 
is also motivated by the composition of its host star. The 
host star 55 Cancri has been reported to be carbon-rich, 
with C/0 = 1.12 ± 0.19 (Mena et al. 2010; E. Del- 
gado Mena, personal communication 2012) and metal- 
rich ([Fc/H] ^ 0.31 ± 0.04; Valenti & Fischer 2005; Mena 
ct al. 2010). However, the stellar C/0 is less than 2- 
sigma above the C/0 of 0.8 required to form abundant 
carbon-rich condensates in protoplanetary environments. 
Furthermore, the C/0 ratio could be overestimated due 
to systematic errors in estimating the C and O abun- 
dances (Fortney 2012). However, if the central value 
of C/0 = 1.12 is taken at face value (see e.g. Carter- 
Bond et al. 2012) and the protoplanetary disk assumed 
to mimic the stellar abundances, the super-Earth 55 Can- 
cri e would indeed be expected to have accreted substan- 
tial amounts of carbon-rich material, such as SiC and C, 
and be depleted in H2O ice. (e.g. Larimer 1975; Bond 
et al. 2008; Johnson et al. 2012). 

In what follows, wc first present, in section [2l inter- 
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Fig. 1. — Mass-radius relations of solid planets with homogeneous composition. The colored curves show mass-radius relations predicted 
by internal structure models of solid planets of uniform composition made of the different species shown in the legend. The black circles 
with error bars show measured masses and radii of known transiting super Earths (Mp = 1 — 10 Mq), adopted from the cxoplanct orbit 
database (Wright et al. 2011). Two values of radii are shown for 55 Cancri e. The red data point shows the radius meas ured in the visible, 
and the blue data point shows a gray radius obtained by combining visible and infrared measurements (see section [2.31 1. 



nal structure models of 55 Cancri e and constraints on 
its interior composition based on its observed mass and 
radius. In section [31 we present thermochemical calcu- 
lations exploring the possible compositions of planetesi- 
mals that may have formed 55 Cancri e. We discuss the 
implications of our results and future work in section |4l 

2. INTERIOR COMPOSITION 

2.1. Internal Structure Model 

We model the internal structure of the super-Earth 
55 Cancri e, and use its observed mass and radius to 
constrain the chemical composition of its interior. Our 
model solves the three canonical equations of internal 
structure for solid planets, namely, (a) mass conserva- 
tion, (b) hydrostatic equilibrium, and (c) the equation of 
state (e.g. Fortney et al. 2007; Seager et al. 2007; Sotin 
et al. 2007; Valencia et al. 2007; Rogers & Seager 2010). 
In each spherical shell of the planet, the equations for 
mass and pressure balance form a system of two coupled 
differential equations, given by: 
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Here, we have adopted the mass interior to each spherical 
shell {Mr) as the independent variable, and we numer- 
ically solve for the radius (r) and pressure (P). The 



density [p) is related to the pressure (P) and tempera- 
ture (T) via the equation of state (EOS) of the material, 
P = P(yO, T). We adopt room-temperature EOSs for the 
species considered in this work, since mass-radius rela- 
tions of solid exoplanets arc known to be insensitive to 
temperature-dependence in the EOS (see e.g. Seager et 
al. 2007). For H2O, MgSiOg (perovskite), Fe, and SiC, 
we use the Birch-Murnaghan (Birch 1952) EOSs com- 
piled in Seager et al. (2007). Our Birch-Murnaghan 
EOS for carbon comprises of the graphite EOS (Wang et 
al. 2012) at low pressures, incorporates the phase transi- 
tion (Bundy et al. 1994) from graphite to diamond at 10 
GPa, and asymptotes to the Thomas- Fermi-Dirac (TFD) 
EOS (Salpeter & Zapolsky 1967) for P > 1000 GPa. For 
diamond, we adopt the EOS from (Dewaele et al. 2008). 
We assume that the planet is spherically symmetric, but 
is vertically differentiated into up to three chemically dis- 
tinct regions, depending on the species under considera- 
tion. The mass fractions of the three regions are given 
as Xi — Mi/Mp] z = [1 — 3], such that xi + X2 -\- X'i = 1. 
Mp is the total mass of the planet and Mi is the mass of 
each of the three layers. 

We solve the system of equations simultaneously using 
a fourth-order Runge-Kutta scheme. The inputs to the 
model are the total mass of the planet [Mp) and the mass 
fractions of two of the three different regions, i.e. xi and 
X2- x^ is obtained from mass conservation, x^ ~ 1 — — 
X2- We start with the upper-boundary {Mr ~ Mp) and 
integrate downward in mass until Mr = 0. The upper 
boundary conditions, which serve as initial conditions. 
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Fig. 2. — Ternary diagrams showing the range of interior compositions allowed by the mass and radii of 55 Cancri e. Two classes of 
interior models were considered, based on the planetesimal compositions predicted by the stellar abundances. Left (a): models composed 
of Fe, SiC, and C, and Right (b): models composed of Fc, MgSiOs, and C. In each case, the red (blue) contours show the constraints 
from the visible (gray) radius (sec Fig. 2). The blue contours are subsets of the red contours. The three axes in each case show the mass 
fractions of the corresponding species. 

sihcate composition would not explain the radius, as 
suggested in previous studies (e.g. Gillon et al. 2012). 
Whereas, if one allows for C-rich compositions, the same 
mass and radius could also be explained by purely solid 
planets with SiC and C based interiors without a volatile 
layer, as can be seen in Figure [TJ 
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and the lower boundary condition is: 



for Mr 
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We use our model in two ways. We first study mass- 
radius relations for planets with homogeneous composi- 
tions. In such cases, xi = 1 and X2 = X3 = 0. We 
choose a grid in Mp, and for each Mp we compute the 
corresponding Rp for the given composition. Since Rp 
is required for the initial conditions but is not known a 
priori, we use a bisection root finder to solve for Rp itera- 
tively such that the lower boundary condition is satisfied. 
Secondly, we also use our model to determine constraints 
on the heterogeneous interior composition of 55 Cancri e 
given its observed Mp and Rp. Here, we explore a fine 
grid in X1-X2 space to constrain the regions that match 
Mp and Rp within their observational uncertainties. 

2.2. Mass-Radius Relations 

Figure [T] shows Mass-Radius relations for homoge- 
neous planets of different compositions (i.e. xi = I, 
X2 = X3 = 0) along with data for several super-Earths. 
Super-Earths with radii above the H2O curve likely host 
gaseous envelopes, since even purely ice compositions will 
not be able to explain their radii. However, those below 
the H2O curve can be explained by various combinations 
of minerals, depending on whether an 0-rich or C-rich 
composition is assumed. An important exception is no- 
ticed for planets such as 55 Cancri e, with radii between 
the silicate and carbon curves. In such cases, assuming 
0-rich compositions would mean that the planet pos- 
sesses a volatile (e.g. H2O) layer, since even a purely 



2.3. Constraints on Interior Composition 

In what follows, we determine the constraints placed 
by the mass and radius of 55 Cancri e on the composition 
space of carbon-rich models of its interior. Interiors of 
carbon-rich planets are expected to be composed of iron 
(Fe), carbon (C), SiC and/or silicates, depending on the 
thermal conditions where the planetesimals originated 
(Kuchner & Seager 2006; Bond et al. 2010; Carter-Bond 
et al. 2012; also see section[3]). Consequently, we consider 
two model families with distinct compositions: Fe-SiC- 
C and Fe-MgSiOa-C. Our planetary interior model in 
each case is comprised of three regions ordered by bulk 
density, an Fe core, overlaid by a layer of SiC or MgSiOs 
(pcrovskite), and an uppermost layer of carbon. In each 
case, Fe-SiC-C or Fe-MgSiOs-C, we constrain the mass- 
fractions of the three regions that explain the observed 
mass and radius of 55 Cancri e. 

Figure [2] shows a ternary diagram of the possible com- 
positions of 55 Cancri e which explain its observed mass 
and radius within their observational uncertainties. We 
adopt the latest estimate of its mass, M = 8.37 ± 0.38 
(Endl et al. 2012), and two values for its radius, 
Rp = 2.04 ± 0.15i?e in the visible (Winn et al. 2011, 
Gillon et al. 2012), and Rp = 2.20 ± 0.12i?e for a gray 
radius obtained by combining the visible and infrared 
(4.5 /.tm) radii (Demory et al. 2011; Gillon et al. 2012). 
The red region (which includes the blue region) in Fig. [5] 
shows constraints from the visible radius, and the blue 
region which is a subset of the red region shows the con- 
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straints from the gray radius. 

Fig. 3a shows the constraints on the mass fractions for 
the Fe-SiC-C class of solutions. A wide range of carbon- 
rich interior compositions can explain the observed mass 
and radius of 55 Cancri e. The visible radius places an 
upper-limit of 41% on the Fe mass fraction, but allows 
a wide range of (SiC, C) combinations. If an Earth- 
like Fe core mass (33%) or less is assumed, the visible 
radius allows a differentiated mantle comprising of all 
possible proportions of SiC and C, constrained only by 
mass conservation (xi + X2 + ~ 1, see section [2.1|) . 
For example, extreme combinations of (Fe, SiC, C) = 
(33%, 67%, 0%) and (33%, 0%, 67%) are both possible. 
However, intermediate compositions such as (33%, 33%, 
34%) may be more representative of the true composition 
since over the long planet formation timescales (Gyrs) 
all the three species are likely to be abundant enough 
to contribute to planetesimal compositions in the inner 
nebular (a < 1.2 AU) where the Fe-SiC-C type of planets 
are more likely to be formed (see e.g. section[3]and Bond 
et al. 2010). 

On the other hand, considering the gray radius nar- 
rows down the solution space, requiring a maximum of 
18% Fe by mass. The SiC and C fractions allowed de- 
pend on the Fe mass fraction. Considering Fe = 10%, for 
example, a minimum of ^70% C is required, with SiC < 
20%. For Fe < 5%, however, all possible combinations of 
SiC and C are allowed, such that SiC + C < 95%. We 
note, however, that the Fc content derived from the gray 
radius may only be a lower limit. The infrared radius 
may contain atmospheric molecular absorption, making 
it slightly larger than the visible radius. Such absorption 
may be caused by an atmosphere (~1000 km) of varied 
compositions, between just ^2 scale heights (Rh) of a 
metal-rich H2-dominated atmosphere and ^30 Rh of a 
CO-dominated atmosphere (see e.g. Fig. 7 of Gillon et al. 
2012). New spectral data would be required to precisely 
constrain this possibility. 

Fig. 3b shows the constraints on the mass fractions 
for the Fe-MgSiOa-C class of solutions. Considering an 
Earth-like Fe content of 33%, the visible radius allows 
for a maximum of 20% MgSiOa, requiring over 47% C, 
which is expected given the dominance of C in the plan- 
etesimal composition in a carbon-rich environment. In 
principle, allowing for lower Fe allows higher MgSiOa. 
For example, a Fe = 20% (10%), allows for MgSiOa up 
to ~60% (80%) and C as low as 20% (10%). Such Fe 
contents are considerably low (compared to Earth's Fe 
fraction of 33 %), and are likely extreme, if not unphysi- 
cal, (also see Gillon et al. 2012), particularly considering 
that the host star has a super-solar metallicity ([Fe/H] = 
0.31 ± 0.04; Valenti & Fischer 2005). However, even a C 
mass fraction of 10% is two orders of magnitude higher 
than that in the Earth which has less than 0.1% C by 
mass (McDonough & Sun 1995). 

3. FORMATION CONDITIONS AND PLANETESIMAL 
COMPOSITIONS 

A carbon-rich 55 Cancri e may also be tentatively con- 
sistent with the composition of the host star. We inves- 
tigate the possible formation conditions of 55 Cancri e 
assuming that the reported stellar abundances represent 
those of its primordial disk (see sections [1] & |4] for chal- 
lenges in this assumption). We adopt the abundances of 



the host star reported by Mena et al. (2010), following 
Carter-Bond et al. (2012). Given the thermodynamic 
conditions in the disk mid-plane, we compute the con- 
densate chemistry in multi-phase chemical equilibrium 
using the commercial software HSC Chemistry. 

We model the protoplanetary disk using a thin disk 
model based on the prescription of Shakura Sz Sunyaev 
(1973), with a turbulent viscosity given by vt = aCg/V,, 
where Cg is the local sound velocity, is the Keplerian 
rotation frequency, and a is the dimensionless viscosity 
parameter. Our model incorporates the opacity law de- 
veloped by Ruden & Pollack (1991) (see e.g. Drouart 
et al. 1999 for details). The temporal evolution of the 
disk temperature, pressure and surface density profiles 
depends upon the evolution of the accretion rate, M, 
which we define (following Makalkin & Dorofeeva 1991) 
to be M = A/o(l + t/to)~^ , where Mq is the initial ac- 
cretion rate, and to is the accretion timescale. We adopt 
S = 1.5, a value that is consistent with the evolution of 
accretion rates in circumstellar disks (Hartmann et al. 
1998). The accretion timescale to is calculated based on 
Makalkin & Dorofeeva (1991) as to = R^j^/ivo, where 
is the turbulent viscosity at the initial radius of the 
sub-disk, Ro- Three parameters constrain Mq and to'- 
the initial mass of the disk Mdq, the coefficient of turbu- 
lent viscosity a, and the radius of the sub-nebula Rd- In 
the present case, we have adopted the parameters used 
by Kavelaars et al. (2011) for their nominal solar nebula 
model, i.e. a = 0.003, i?Do = 15 AU, and Mdo = 0.06 

Mq. 

The refractory condensate composition as a function of 
the temperature (T) is shown in Fig. [3l At high temper- 
atures (T > 1100 K), SiC and Fe are the most dominant 
condensates. As the temperatures cool to between 700 
and 1100 K, pure carbon (as graphite) becomes the most 
dominant condensate, followed by Fe and silicates, e.g. 
Mg2Si04. However, at low temperatures, T below 700 K, 
Mg2Si04 dominates the condensate composition despite 
the high C/0 ratio (also see Bond et al. 2010; Carter- 
Bond et al. 2012). Figure 1 also shows representative 
orbital separations (a) in the proto-planetary disk and 
their corresponding temperatures at different epochs in 
the disks evolution. For a < 1.2 AU in the disk, SiC, 
Fe, and C condense early in the disks lifetime, and hence 
dominate the planetesimal composition. For a ~ 1.2 - 
1.6 AU, graphite far dominates the planetesimal compo- 
sition, followed by Fe and Mg2Si04. At lower tempera- 
tures, however, for a > 1.6 AU, Mg2Si04 dominates the 
planetesimal composition. 

Depending on whether the planetesimals that formed 
55 Cancri e originated primarily within or beyond ^1.2 
AU, its interior can bear end-member compositions of Fe- 
SiC-C or Fe-silicates-C, respectively, or an intermediate 
composition. Since the protoplanet likely formed farther 
than the planet's current orbital separation (0.015 AU) 
and migrated in, it likely accreted substantial amounts of 
Fe-rich and C-rich condensates irrespective of its original 
location. On the other hand, it is unlikely that the planet 
formed in-situ because the disk temperatures at this lo- 
cation prevent the condensation of dominant refractory 
condensates. 

4. DISCUSSION 
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Fig. 3. — Chemical composition of refractory condensates expected in the proto-planetary disk of the 55 Cancri system. The colored 
curves show molar mixing ratios of the major species, shown in the legend, in thermo-chemical equilibrium as a function of temperature in 
the proto-planetary disk-midplane. A representative pressure of 10~ bar is assumed. The elemental abundances of the host star (Mona 
et al. 2010; Carter-Bond et al. 2012) are assumed as initial conditions, and the equilibrium computations were performed using the HSC 
Chemistry software. The black filled circles at the top of the figure show the temperatures of the disk at different orbital separations and 
ages. 



We show that the sum-total of available data, includ- 
ing the planetary mass and radius, and the stellar abun- 
dances, support the possibility of a carbon-rich interior 
in 55 Cancri e. Numerous studies have suggested the pos- 
sibility of carbon-rich cxoplanets, both theoretically (e.g. 
Kuchner & Seager 2006, Bond et al. 2010, Madhusud- 
han et al. 2011b) and observationally (Madhusudhan et 
al. 2011a, Madhusudhan 2012). However, the prevalent 
practice is to assume Earth-centric compositions, com- 
prising of Fe, silicates, and H2O, in explaining super- 
Earth observations. A carbon-rich 55 Cancri e would 
represent a departure from Earth-centric mineralogies in 
rocky exoplanets. Nevertheless, follow-up efforts are re- 
quired to confirm a carbon-rich 55 Cancri e. 

High confidence elemental abundances of the host star 
can help refine the constraints on the formation condi- 
tions of 55 Cancri e. Given currently reported abun- 
dances, a C/O>0.8, required to form C-rich refractories 
(Larimer 1975), holds only at ~2-tT and the C/0 ra- 
tio may also have been overestimated (Fortney 2012). 
Furthermore, we assume that the protoplanetary disk 
bears the same composition as the host star. While this 
may be a reasonable assumption for deriving refractory 
compositions of planetesimals forming rocky planets (see 
e.g. Bond et al. 2010; Carter-Bond et al. 2012; John- 
son et al. 2012), the possibility of local perturbations 
in the C/0 ratio of the disk may not be entirely ruled 
out; i.e. the possibility that carbon-rich planets may 
form around oxygen-rich stars and vice versa (Kuchner 
& Seager 2006; Madhusudhan et al. 2011b; Oberg et al. 



2011). In particular, carbon-rich gas giant planets may 
form by predominantly accreting effectively carbon-rich 
gas in regions beyond the water-ice condensation line in 
an oxygen-rich disk (Oberg et al. 2011). 

High-precision spectroscopy in the future may also 
be able to constrain the presence of an atmosphere in 
55 Cancri e and its composition. Currently, both re- 
ported radii, visible and infrared, are consistent at the 1- 
tr uncertainties, though the latter is slightly larger (Gillon 
et al. 2012). Excess infrared absorption may be caused 
by molecular species and/or collision-induced absorption 
(Miller- Ricci et al. 2009) which could be constrained by 
future observations. 

A carbon-rich interior in 55 Cancri e has far-reaching 
implications for geophysical processes in its interior. The 
Earth's interior is oxygen-rich, with carbon less than 
0.1% by mass (McDonough & Sun 1995), in contrast to 
a carbon-rich interior of 55 Cancri e. Both carbon and 
SiC have over two orders of magnitude higher thermal 
conductivities than the silicates and oxides that dom- 
inate the Earth's mantle (Goldberg et al. 2001), im- 
plying that the thermal evolution of 55 Cancri e may 
be different from that of oxygen-rich planets. The high 
thermal conductivities, compounded with the paucity of 
water, may influence the likelihood of plate tectonics 
in carbon-rich super-Earths (Valencia et al. 2007; Ko- 
renaga 2010), though new measurements of viscosity of 
carbon-rich materials at high pressures and high tem- 
peratures are required to constrain this possibility. A 
carbon-rich 55 Cancri e bolsters the idea that rocky ex- 



6 



oplancts may be diverse in their chemical compositions 
and, hence, in their geophysical conditions. 
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